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We discuss theoretically, using density-functional theory, the phase stability of nematic and smec- 
tic ordering in a suspension of platelets of the same thickness but with a high polydispersity in 
diameter, and study the influence of polydispersity on this stability. The platelets are assumed to 
interact like hard objects, but additional soft attractive and repulsive interactions, meant to repre- 
sent the effect of depletion interactions due to the addition of non-absorbing polymer, or of screened 
Coulomb interactions between charged platelets in an aqueous solvent, respectively, are also consid- 
ered. The aspect (diameter to thickness) ratio is taken to be very high, in order to model solutions 
of mineral platelets recently explored experimentally. In this regime a high degree of orientational 
ordering occurs; therefore the model platelets can be taken as completely parallel and are amenable 
to analysis via a fundamental-measure theory. Our focus is on the nematic vs. smectic phase in- 
terplay, since a high degree of polydispersity in diameter suppresses the formation of the columnar 
phase. When interactions are purely hard, the theory predicts a continuous nematic-to-smectic 
transition, regardless of the degree of diameter polydispersity. However, polydispersity enhances 
the stability of the smectic phase against the nematic phase. Predictions for the case where an 
additional soft interaction is added are obtained using mean-field perturbation theory. In the case 
of the one-component fluid, the transition remains continuous for repulsive forces, and the smectic 
phase becomes more stable as the range of the interaction is decreased. The opposite behaviour 
with respect to the range is observed for attractive forces, and in fact the transition becomes of first 
order below a tricritical point. Also, for attractive interactions, nematic demixing appears, with 
an associated critical point. When platelet polydispersity is introduced the tricritical temperature 
shifts to very high values. 

PACS numbers: 61.30.Cz, 61.30.Hn, 61.20.Gy 



I. INTRODUCTION 

The issue of polydispersity is crucial to understand 
phase behaviour in experiments on colloidal suspensions 
of spherical and anisometric particles. The latter, of ei- 
ther rod- or plate-like shape, are known to form liquid- 
crystalline phases and their phase equilibria is 
largely affected by size dispersity |9l-[ll| . Colloidal parti- 
cles can never be made truly identical, and even a small 
amount of polydispersity is unavoidable in real samples. 

Platelets with varying degrees of polydispersity, in par- 
ticular, are being extensively studied in recent years. In 
this material, polydispersity greatly facilitates gelation, 
and in fact the nematic phase is difficult to reach in sus- 
pensions of charge-stabilised platelike particles (the most 
common plate-like particles) at particle concentrations 
below that at which a glassy state is formed How- 
ever, it was possible to obtain a nematic phase in hard- 
like platelets, e.g. gibbsitc particles with steric stabilis- 
ers [y] . It has been demonstrated that the observation of 
equilibrium nematic phases in suspensions of charged col- 
loidal platelets requires fine-tuning of the ionic strength 
of the solvent • Therefore, the concomitant effects 



of polydispersity and soft interactions seem to be an im- 
portant issue as regards the stabilisation of the nematic 
phase. 

Besides the nematic phase, liquid-crystalline phases 
with partial spatial order in colloidal suspensions of 
platelets are now being investigated 1J| . The most com- 
mon phase is the columnar phase, which admits a sur- 
prisingly large (up to 25%) degree of polydispersity in 
diameter. Suspensions of sterically-stabilised gibbsite 
platelets with added non- absorbing polymer exhibit a 
rich phase diagram, with gravity-induced three-phase co- 
existence regions involving isotropic, nematic and colum- 
nar phases, and enhanced fractionation effects By 
contrast, observations of the layered smectic phase are 
very rare in these systems 16-1^, since the columnar 
phase is very stable at high particle concentrations. 

In a recent paper [l^, smectic ordering was observed 
in colloidal suspensions of equally-thick charged platelets 
with high polydispersity in diameter and in a solvent with 
low ionic strength. It is though that the extreme di- 
ameter polydispersity (> 30%), possibly combined with 
the monodispersity in thickness, is responsible for the in- 
stability of the columnar phase in favour of the smectic 
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phase. It would be interesting to theoretically analyse the 
effect of polydispersity on the formation of the spatially- 
ordered smectic phase in systems where the columnar 
phase is suppressed. 

On the other hand, additional soft interactions are 
known to affect the phase diagram. These interactions 
can be screened Coulomb particle interactions due to 
the addition of salt to an aqueous solution of charged 
particles, or attractive depletion forces [l^ between the 
colloidal particles arising from addition of non-absorbing 
polymer to a suspension of effectively hard platelets [20l - 
|23| . The possibility of theoretically predicting how both 
factors (polydipersity and soft interactions) influence the 
formation of phases with partial spatial order (smectic 
and columnar phases) would be highly desirable. 

In the present paper we study the effects of polydis- 
persity and soft interactions on the nematic vs. smectic 
phase equilibrium. The work was inspired by the pa- 
per by Sun and coworkers [3], who synthesised novel 
platelet particles from a-ZrP minerals by a special exfo- 
liation procedure, which creates a perfectly monodisperse 
thickness distribution [24,, ,25;] with high polydispersity in 
diameter. Low ionic-strength aqueous solutions of these 
charged particles form equilibrium nematic phases, fol- 
lowed by an incipient smectic ordering which equilibrates 
very slowly; no hint of columnar ordering was found (isj . 

In the experiments by Sun et al. [Tsj the ionic strength 
of the solution was not controlled. However, addition of 
salt to the suspension will modify the long-ranged repul- 
sive interactions between platelets. Our second aim is to 
predict how the nematic vs. smectic phase equilibrium of 
the suspension will be modified when the range of inter- 
actions is changed. These results could also apply to sus- 
pensions of colloidal platelets with non-absorbing poly- 
mer, where effectively attractive depletion force between 
platelets usually result in demixing between phases with 
different concentrations of colloids and polymers [26| : 
therefore, we also consider attractive soft interactions of 
varying range. As shown in Ref. [27| . addition of non- 
absorbing polymer to a one-component fluid of hard rods 
may give rise to isotropic liquid-liquid phase separation 
ending in a critical point (for large polymer coils) , and to 
the presence of a broad coexistence region of isotropic- 
nematic, isotropic-smectic and isotropic-crystal phases. 
Binary mixtures of colloidal disks and polymers exhibit 
not only isotropic-isotropic demixing, but also nematic- 
nematic demixing 28] . However nonuniform phases were 
not taken into account in the latter study. As will be 
shown below, soft interactions between highly oriented 
platelets can induce nematic-nematic demixing and a 
dramatic broadening of the nematic-smectic transition 
region for sufficiently long-ranged interactions. 

The effect of polydispersity on the phase behavior of 
hard rods has been theoretically studied at the level 
of density functional theory (using the Onsager second 
virial theory) mainly for uniform liquid crystal phases 
(isotropic or nematic) (29l - [3l| . Also, a bimodal polydis- 
perse distribution function in aspect ratio for hard board- 



like particles was used to study the effect of polydisper- 
sity on the phase stability of biaxial nematic vs. nematic- 
nematic demixing [s^, [l^l • The limit of zero polydisper- 
sity gives rise to a perfect binary (bidisperse) mixture, a 
case studied by Varga and coworkers jM] using Onsager 
theory. Harnau et al. [3^ and Bier et al. [36| have also 
studied perfect binary mixtures of platelets using density- 
functional theory [361 ] in the restricted-orientation ap- 
proximation. Finally, a simulation study was done by 
Bates and Frenkel on the influence of polydispersity on 
nonuniform phases of rod-like particles: it was shown 
that the smectic phase of hard colloidal rods can be 
destabilized with respect to the columnar phase for high 
enough length polydispersity. 

Density-functional studies of the relative stability of 
nonuniform (smectic, columnar and crystal) phases in 
polydisperse fluids made of anisotropic particles consti- 
tute a great challenge. The main difficulty is related 
to the numerical minimization of the functional as the 
density profile is a function of the spatial, angular and 
polydisperse variables. The present work is a first at- 
tempt to tackle the study of smectic phases in diameter- 
polydisperse platelets using density-functional tools and, 
in particular, fundamental-measure theory (FMT) of par- 
allel hard cylinders [37| . 

We first give a brief overview of the theory in Sec. [Hi 
In Sec. mil we examine the influence of polydispersity on 
the phase behaviour of hard parallel platelets, while the 
case of soft platelets is treated in Sec. IIVI We end with 
some conclusions in Sec. El 



II. MODEL AND THEORY 

The particle model consists of cylinders with paral- 
lel axes of revolution. This approximation is justified in 
the regime of high packing fractions where the nematic- 
smectic (NS) transition occurs (obviously the isotropic- 
nematic transition cannot be treated within this scheme, 
but this is not the aim of the present study). Cylinders 
have the same thickness, Lq, in line with the particles 
obtained by Sun et al. [la] , but are polydisperse in diam- 
eter. We assume the fluid to consist of a multicomponent 
mixture of species with different radii R. In the thermo- 
dynamic limit the fluid contains infinitely many species, 
and polydispersity can be characterised by means of a 
continuous size distribution. The particle radius will be 
characterised by the scaled variable r = R/Rq, with Rq 
the mean radius. The number fraction of a given species 
with radius r will be p{r), the (continuous) radius distri- 
bution to be specified later. 



A. Density functional theory 

The statistical mechanics of a fluid of N parallel cylin- 
drical platelets in a volume V will be investigated using 
the FMT density-functional approach for a multicompo- 
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FIG. 1: Schematic of the parallel cylindrical particles used in 
this work. 



ncnt mixture of parallel cylinders presented in [37| . Here 
we extend this approach to the polydisperse case. 

We choose the cylinder axes, the nematic director and 
the layer normal in the smectic phase to be parallel and 
along the z direction. In the density-functional approach 
the central quantity is the number density p(z, r), which 
gives the local density at z of the species with radius r 
and satisfies 



dr I drp{z, r) = N. 
'v 



(1) 



Note that r — {x, y, z) refers to the space coordinate 
vector while r is the scaled radius of particle (not the 
absolute value of r). The excess part of the local free- 
energy density (assuming smectic symmetry) is 



<i>(z) = PfUz) 

= -710 (1 - "3) + 



1 — 713 



r, (2) 



where is the local excess free-energy density, (3 — 
1/kT, and k is Bolztmann constant, with T the temper- 
ature. The functions no(z), n\{z), ni{z), n\{z), n^(z) 
and 713(2) [in ^ their z dependence has not been ex- 
plicitly indicated] are obtained from the number density 
p{z,r) using the expressions: 



no{z) 



Xo*'^^°^l {z), Ti\{z)^ t:\Mq *uj'~^'^\ {z), 



nt{z) = Mi*c^(3) nl{z)^ M2 * w^"^ (z). 



,(0) 



ni{z) 



(z), 773(2) = 7r 



A^2 * w 



(3) 



(3) 



where the symbol * stands for convohition. The functions 
Xa(z) = / drp{z,r)R'', (4) 



are generalized moments of the polydispersity distribu- 
tion, while 



/3)(z) = e(|^- 



,(5) 



with S and Q the Dirac delta and Heaviside functions, 
respectively. The excess free-energy functional is ob- 
tained by integration over the system volume V, i.e. 



/3J-"ox[p] = Jydr^{z). Adding the ideal contribution 
PTidlp], with 



dr / drp(z, r) {log [p(z, r)Ar] - 1} , (6) 
Jv 



where is the thermal wavelength of the species with 
radius r, we obtain the total free-energy functional as 
J^[p] — J-id[p] + -^ox[p]- The equilibrium state of the sys- 
tem follows from functional minimisation of J-[p\ with 
respect to the local number densities p{z,r). 

The radius distribution function p{r) was chosen to be 
a Schultz distribution. 



P{r) = exp [-(7/ + l)r] 



(7) 



where 1/ is a free parameter that controls the width of 
the distribution. Once p(r) is chosen, the local density 
p{z,r) must satisfy the normalisation relation 



dvp{z, r) 



dr / drp{z,r) 
V Jq 



^J^drp{z,r). 



(8) 



Defining the moments as 



drp{r)r° 



(9) 



the distribution (|5]) fulfills uiq — mi — 1. The polydis- 
persity coefficient A, defined as usual as the width of the 
radius distribution, A = (7712 — 771^)^/^/7711, is related to v 
by A = (t^-I- 1)~^/^. In the smectic phase, local moments 
ma{z) of the density distribution can also be defined: 



1 /"^ 

ma{z) = — / drr'^p{z,r), 
Po Jq 

where po = N/V is the mean density. Obviously 
If 1 

— / dvniaiz) = — dznialz) = 771^, 
V Jv d Jq 



(10) 



(11) 



IV " Jo 

where d is the smectic period. 



III. HARD PLATELETS POLYDISPERSE IN 
DIAMETER 

In the one-component case (platelets of the same di- 
ameter) the NS transition is located at a value of pack- 



ing fraction rj. 



(0) 



0.314 [38|, where the packing frac- 



tion is defined as 77 = p^v, with v — ttR^Lq the mean 
particle volume and po = N/V the mean density. The 
transition is continuous. The question we would like to 
answer is whether the transition is still continuous when 
polydispersity is introduced and how the packing fraction 
at the transition behaves as the sample polydispersity is 
changed. 
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FIG. 2: Packing fraction rjs at the NS spinodal, relative to 
that of the monodisperse fluid rji''\ as a function of polydis- 
persity parameter A for hard platelets. Labels indicate the 
nematic (N) and smectic (S) regions. The inset is the radius 
distribution function at the spinodal for A = 0.6 




A. Nematic-smectic bifurcation 

As a preliminary step prior to the full minimization of 
the functional, we have calculated the spinodal (or bi- 
furcation) line for the NS transition of the fluid. This 
can give us an indication of the basic trends. The cal- 
culation starts by perturbing the nematic with a small- 
amplitude density wave, p{z,r) = pop(r) + e(r)cosfcz, 
where k — 2n/d is a wavenumber. The fluid response to 
the perturbation is obtained from the Fourier transform 
of the direct correlation function (taken from the sec- 
ond functional derivative of the excess free-energy func- 
tional): 

-poc(fc;r,r') = co(fc)(r + r')^ 

+ ci(fc)rr'(r2 + r'^) + C2{k){rr')^{12) 

where r, r' are the scaled radii of two platelets. Here we 
have defined the coefficients 

coik) = y[2x{k) + ym2x\k/2)] , 

ci{k) = 2y2[2x(fc) + (l + 2yTO2)x'(fc/2)], 

C2(fc) = y'[2{l + 2y)x{k) 

+ (1 + 27/(2 + ma) + 6TO22/')x'(fc/2)] , (13) 

where x(fc) = sinfc*/fc* (with k* — kL^ the scaled wave 
number) and y = rj/ (1— r/ma). The instability of the fluid 
against the density wave is obtained from the eigenvalue 
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FIG. 3: (a) Scaled pressure pPv, and (b) scaled free energy 
j3Fv/V as a function of packing fraction rj for the nematic 
(dashed curve) and smectic (continuous curve) branches of 
hard polydisperse platelets. The circle indicates the location 
of the bifurcation point. The value of polydispersity is A = 
0.447. 



problem 



i-{r)=pop{r) I dr'c{k,r,r')e{r'), (14) 



from which the wavenumber k^ and packing fraction r]s 
at bifurcation can be obtained. 

Fig. [5] shows the packing fraction of the NS spinodal, 
?7s, relative to that of the monodisperse fluid, ri'f \ as a 
function of polydispersity A. The inset shows the radius 
distribution function for the particular value A = 0.6. 
The packing fraction at bifurcation decreases as particles 
become more polydisperse in radius, so the smectic phase 
becomes more stable with respect to the nematic with 
increasing polydispersity. This trend is opposite to that 
observed in fluids of hard rods, where polydispersity in 
length postpones the onset of smectic stability due to 
incommensurability of the particle lengths with a smectic 
period. 
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FIG. 4: Profiles of tlie moments ma{z), witli a — 0,1,2 
(continuous, daslied and dotted curves, respectively), along 
one smectic period for hard polydisperse platelets in the case 
A — 0.333 and -q — 0.455. The smectic period is d = 1.212Lo. 



B. Full minimisation and search for coexistence 

A full minimisation of the functional is required in or- 
der to gain more insight into the true nature (i.e. contin- 
uous vs. first-order) of the NS transition and to obtain 
information on the equation of state, smectic period, and 
possible fractionation effects in the case A > (i.e. the 
fact that the radius distribution function can be different 
in the two coexisting phases). 

We wish to calculate the density profile p{z,r) of a 
stable smectic phase with a fixed number density po and 
a prescribed radius distribution function p{r). Note that 
the nematic phase will be unstable or metastable. In 
the minimisation of the free-energy functional per unit 
volume, I3J-/V, the density distribution is subject to the 
constraint pop{r) = {p{z,r))j^ [see Eqn. with the 

shorthand notation(- • ■ )^ = Jq{- ■ ■) for an average 
over the smectic period. Direct functional minimisation 
leads to 



p{z,r) 



PoPir) exp [c(^)(z,r)] 
(exp[cW(z,r)])^ 



(15) 



with c^^\z,r) the one-body correlation function. This 
function is a quadratic polynomial in r, with coeffi- 
cients which are functionals of the moments ma{z) = 
p^^ drr°'p{z,r). From pS]) we obtain a set of self- 
consistent non- linear integral equations in ma{z): 



ma{z) = 



dr 



r"p[r) exp [c^"'"' (z, r)] 
(exp [c'-^Hz,r)])^ 



, a = 0,l,2.(16) 



Expanding the local moments ma{z) in a Fourier series 
containing the wavenumber k = 271 /d^ Eqns. (jl6p can be 
written as self-consistent equations for the corresponding 
expansion coefficients. The smectic period d is obtained 



from the free energy by minimisation. For high pack- 
ing fraction this process gives moments •ma{z) that are 
not constants: this corresponds to the smectic branch. 
Otherwise the nematic branch is obtained. 

Next we briefly outline how the conditions for NS co- 
existence were obtained. In equilibrium the chemical po- 
tentials of each species should be equal in both phases: 
^(*)(r) = /i(°)(r) {i = N,S for nematic and smectic, re- 
spectively). /i'^°-'(r) is calculated from the level-rule con- 
straint 



lP^{r) + (1 - l)Ps{r) = Popir), 



(17) 



which implies the conservation of the total number of 
particles; po and p{r) are the number density and radius 
distribution function of the parent phase. < 7 < 1 
is the fraction of the total volume occupied by the N 
phase. Note that the number density of the S phase 
is Ps{r) = {ps{z,r)) ^. Using the definition for chemi- 
cal potential, p'^^\r) = 5T / 5pi{z,r)\ei^, the equilibrium 
conditions together with the level rule provide a set of 
equations for the moments mi"'' and rna\z): 



,(n) 



drr^pir) exp [cl,^^ (r)j T{r), (18) 



m^"^(z)= / drr"p{r)exp ci^\z,r) T{r), (19) 
'0 



T~"^(r)=7exp c';!'\r) 

+ (l-7)(exp[4^)(z,r)])^. (20) 

These equations, together with the condition for mechan- 
ical equilibrium, i.e. equality of pressures of both phases 
p(s) ^ p(N) jij^ Q^j. ^q^qI pp = {d<t>/dn3)^], provide 

equilibrium moments of each phase and the number den- 
sity Po of the parent uniform phase at coexistence. 

A typical feature of polydisperse systems is that the 
composition of the two coexisting phases varies as the 
phase transition takes place. In our case this means the 
following. Let us fix the value of A (the so-called dilution 
line) . Then, as the total packing fraction rj of the parent 
nematic phase is increased, first an infinitesimal amount 
of smectic material, called shadow smectic, will appear, 
coexisting with the cloud (parent) nematic phase. As ry 
is further increased, the amount of smectic material will 
grow, and eventually the nematic will disappear. The 
opposite process occurs by starting from a parent smectic 
phase and decreasing 77; a smectic cloud point will be 
reached when the first (shadow) nematic material can be 
observed in the sample. 

The cloud-nematic and shadow-smectic coexistence 
curves (which are more relevant experimentally since one 
usually starts from a parent nematic phase), can be cal- 
culated from Eqs. ^ and ^ by setting 7=1. This 
results in 



drr°'p(r)e^ 



a = 0,1, 2, (21) 



6 



where Ac(^^(z, r) = ci^\z, r) — ^(r) is the difference in 
the one-body correlation functions of the two phases. In 
this case the nematic radius distribution function coin- 
cides with p{r) (that of the parent phase) . The smectic- 
cloud and nematic-shadow calculations can be imple- 
mented by setting 7 = in Eqs. (fTS]) and (|20p . In what 
follows, all coexistence results correspond to the cloud- 
nematic and shadow-smectic curves for different polydis- 
persities. 

Our first result of the full minimisation is shown in 
Fig. [31 where the scaled pressure (3Pv and the scaled 
free energy pPv/V are plotted as a function of packing 
fraction rj for the case A = 0.447. Calculations using 
many different initial guesses to solve the self-consistent 
equations always lead to a single stable solution, which 
means that the nematic branch always bifurcates tangen- 
tially to a smectic branch. This implies that there is no 
coexistence, and that the transition is always continuous. 
The scenario is the same for polydispersities larger than 
the one used in the calculations shown in the figure. 

It is instructive to examine the structure of the smectic 
phase by looking at the profiles of the moments ma{z). 
Fig. 21 plots the moments in one smectic period for the 
case A = 0.333 and packing fraction rj = 0.455, i.e. well 
inside the smectic region. As can be seen from the figure, 
the zeroth moment, mo(z), reflects the structure of the 
total density, i.e. pronounced peaks located at the smec- 
tic layers with a period of d. The first moment, mi(z), 
gives information about the mean particle radius at each 
location; at the layers the mean radius is ~ 1.65i?0j a-nd 
it decreases away from the layers: larger particles are 
located preferentially at the layers, whereas smaller par- 
ticles tend to stay at the interstitials. 

The latter effect can be studied in more detail by 
analysing the density profiles p{z, r), along with the local 
radius distributions 



r) 



p{z,r) 



drp{z, r) 



(22) 



as a function of r for given z. This is shown in Figs. El^a) 
and (b) for four different values of z. Complementary 
plots are presented in Figs. WLc) and (d), where the same 
functions are given as a function of z for three particular 
values of r. The conclusion that can be drawn from the 
figures is that there exists a microsegregation effect in the 
smectic phase. Specifically, particles with larger radii 
are preferentially located at the smectic layers, i.e. at 
z = or d, and possess a wider distribution function; 
by contrast, at the interstitial, z = d/2, the distribution 
function h{z,r) has a maximum at smaller values of r 
and is much more narrow. These microsegregation effects 
become more pronounced as the degree of polydispersity 
is increased. 

The increased stability of the smectic phase in fiuids 
of hard polydisperse platelets, shown in Fig. [21 is due 
to the larger packing efficiency (lower excluded volume) 
that platelets can achieve at the smectic layers when 
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FIG. 5: Scaled density distribution p{z, r)v and normalised 
density profile h{z, r) for a smectic phase of polydisperse hard 
platelets with packing fraction r; = 0.452 and A = 0.523. In 
(a) and (b), distributions are plotted as a function of the 
scaled radius r for particular values of the z coordinate, given 
as labels (for reference, the smectic period is d = 1.211I/o). 
In (c) and (d), distributions are plotted as a function of the 
z coordinate for particular values of r. 



there is a continuous distribution of radii. As the dis- 
tribution becomes wider, and through a microsegrega- 
tion mechanism, larger packing fractions can be achieved 
at the quasi 2D-layers, which leads to a stabilisation of 
the layered arrangements with respect to the uniformly 
distributed platelet configurations in the nematic phase. 
In fluids of polydisperse rods, microsegregation also ex- 
plains the destabilisation of the smectic phase, since par- 
ticles of different lengths cannot be arranged favourably 
into identical layers because of inefficient packing, and 
the minority species is expelled to the interstitials. 



IV. SOFT PLATELETS 

Next we consider hard platelets that interact via an ad- 
ditional soft potential. The soft potential will be treated 
by means of the usual mean-field approximation, neglect- 
ing correlations. First we treat the case of a monodis- 
perse fluid, leaving the general case of a polydisperse fluid 
for the last section. 



A. Monodisperse fluid 

The soft potential should reflect the anisotropics of 
particle interactions. We have chosen a functional form 
corresponding to a modifled Yukawa potential where 



7 



transverse and longitudinal coordinates are scaled with the particle sizes R and Lq, respectively: 
I 



Vsoft{r±,z) = =Fe < 



[ 0, 



v/(r±/2i?)2 + (z/Lo)2-l 



r 



\z\ > Lo, 
or 

\z\ < Lq and r± > 2R, 
otherwise, 



(23) 



where e > 0. Here r± is the interparticle relative distance 
(measured from the platelet centres) in the xy plane, z 
the relative distance along z, and A an inverse interaction 
range parameter. Scaling the coordinates as in (j23p is a 
natural choice, as it produces ellipsoidal-like equipoten- 
tial surfaces and discriminates between face-to-face and 
side-by-side configuration of two platelets, which should 
have very different interparticle forces. 

The effect of the soft interaction on the free-energy 
functional is obtained from perturbation theory in the 
mean-field approximation: 



- 2d / 



dz'p{z)p{z')u{\z' - z\) (24) 



with the effective potential 

u{z) = 27r / dr^r^/314oft {r± , z) . (25) 
Jo 

We use a Gaussian parametrization for the smectic den- 
sity profile: 



P(z) = pd(^) ^ ^ exp [-a(z - nd)2] , 



(26) 



where n is the set of integer values: 0, ±1, ±2, ... 

In the present case of monodisperse particles, the pro- 
cedure outlined in Sec. IIII BI for obtaining the equilibrium 
structure of the fluid can be made much simpler. Thus, 
for given values of scaled temperature T* = fcT/e and 
interaction range parameter A, the calculation proceeds 
by fixing the value of the packing fraction 77, and then 
minimising the free-energy with respect to the Gaussian 
width a and the smectic period d to find the equilibrium 
density profile. Having obtained the free-energy branches 
as a function of r/ for both nematic and smectic phases, 
the double tangent construction is used to find the coex- 
isting packing fractions. Changing T and repeating the 
whole process, the phase diagram of the model in the 
T — rj plane is obtained. Considering the cases — or -I- in 
(|23p one can cover both attractive and repulsive interac- 
tions. 

We now turn to the results. In the case of attractive 
forces, the NS transition can be of first order or contin- 
uous, depending on the temperature. This is shown in 



Fig. ini for two values of the inverse range parameter A. 
There is a tricritical point (indicated by an open circle in 
the figure), below which there is NS coexistence. Above 
the tricritical temperature the transition is continuous. 
This scenario is independent of the value of the interac- 
tion range. However, an interesting feature of the phase 
behaviour is that, when the interactions are sufficiently 
long-ranged (i.e. A sufficiently small), the nematic free- 
energy branch presents a region of instability, indicat- 
ing the existence of nematic-nematic (Ni-fN2) demixing 
below a critical temperature [with an associated critical 
point indicated by a shaded circle in Fig. [DJa)]. Obvi- 
ously, as the inverse length decreases, nematic-nematic 
demixing becomes even more pronounced. This is shown 
in Fig. [Tl^a), where the N1-I-N2 demixing region in the 
T — rj plane is represented for different values of the in- 
verse lengths, A = 0.5, 1 and 2. The effect of the range 
of the interactions on the smectic stability, which is rel- 
atively weak, can be appreciated from the change in the 
location of the NS spinodal as A is varied. As can be seen 
from Fig. [Tja), smectic ordering is slightly favoured as 
the range of the interactions is increased. 

The case of repulsive interactions is quite different. No 
nematic demixing transition exists in this case of repul- 
sive interactions [Fig. Efb)], and the NS transition is al- 
ways continuous. Again the effect of varying the inverse 
length is weak, but different from the attractive case: for 
attractive forces, smectic ordering is favoured as the in- 
teraction range is increased, but in the case of repulsive 
forces, longer-ranged interactions disfavour the formation 
of a layered structure with respect to a spatially uniform 
configuration. In any case, the effect of the interaction 
range on the NS transition is larger for repulsive forces. 



B. Polydisperse fluid 

When platelets are not equal in size several modifica- 
tions in the theory are necessary. The first is that density 
profiles depend on the polydispersity variable r, and the 
soft potential Vso{t{fL, z, r, r') also contains a dependence 
via the radii of the two interacting platelets. The mean- 
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FIG. 6: Phase diagram of the monodisperse attractive 
Yukawa fluid. Continuous curves: coexistence boundaries for 
the nematic-smectic or nematic-nematic transitions. Dashed 
curves: nematic-smectic spinodal. Open circle: critical 
points. Shaded circle: tricritical point. Dotted horizontal 
line: tieline of triple-point N1-I-N2-I-S coexistence, (a) A = 1. 
(b) A = 2. 



J 



/3-^mf[p] 



nd /'oo 

Yd ' 



dz' I dr dr' 
Jo 



X p{z,r)p{z',r')ui\z-z'\,r,r'), (27) 



where u{z, r, r') is defined as in (|25l) . Next the soft inter- 
action potential has to be specified. The soft potential 
will be defined as an extension of the one used previ- 
ously, Eqn. (|23p . but for two platelets of different radii 
R = rRo, R' = r'Ro: 



Vsoitirj_,z,r,r') = -e < 



r^/{R + R')f + {z/Lof-l 



r^/{R + R')f + {z/Lof 



[0, 



\z\ > Lq 
or 

\z\ < Lq and r± > R + R' 
otherwise 



(28) 



Again, since we are now considering a polydisperse fluid, 
a proper calculation of the cloud nematic and shadow 
smectic curves has to be implemented. The radius distri- 
bution function used in the following contains a Gaussian 
tail: 



(29) 



2s-:. 



r[(i. + i)/2]' 



r[(^ + 2)/2] 
r[(^. + i)/2]- 



(30) 



This choice is related to the better convergence reached 
in the numerical implementation of Eq. (|2ip . 

We have performed coexistence calculations and ob- 
tained complete phase diagrams including nematic and 
smectic regions of stability, in this case only for the at- 
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FIG. 7: Scaled temperature-packing fraction T* — rj phase 
diagrams of the monodisperse model with (a) attractive, and 
(b) repulsive soft interactions. Nematic (N) and smectic (S) 
phases are labelled by symbols. In the case of demixing, two 
nematic phases, Ni and N2, coexist below a critical point, 
indicated by circles. In both panels the inverse interaction 
length is A = 0.5 (dashed curves), 1.0 (continuous curves) 
and 2.0 (dotted curves). 
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FIG. 8: Phase diagram of the polydisperse attractive Yukawa 
fluid. Continuous lines: coexistence boundaries for the poly- 
disperse fluid with A — 0.294. Dashed lines: monodisperse 
fluid. Dotted line: nematic-smectic spinodal. (a) A = 1. (b) 
A = 2. 



tractive Yukavifa-like potential (polydispersity is not ex- 
pected to affect the phase behaviour significantly in the 
case of repulsive forces). The phase diagrams are shown 
in Figs. Ufa) and (b), the first for an inverse length A = 1 
and the second for A = 2. For the sake of comparison, 
the results for the monodisperse fluid have been superim- 
posed. As far as N1-I-N2 coexistence is concerned, both 
the cloud and shadow curves of each nematic phase were 
calculated; however, their difference cannot be appreci- 
ated at the scale of the figure, so only the cloud curve is 
plotted (this curve chosen for consistency, as in the case of 
the NS coexistence only the cloud nematic-shadow smec- 
tic curve was computed). The most important change 
brought about by polydispersity is that the tricritical 
point moves to very high values even for small values of 
A, and therefore the NS phase transition becomes of first 
order in the reasonable range of temperatures. Again, 
for the longer-ranged potential, there exists an region of 
nematic demixing that is missing for shorter-ranged in- 
teractions, a feature not much affected by polydispersity. 



To show the effect of polydispersity on the NS transi- 
tion in attractive platelets and ascertain how the nature 
of the transition changes from the mono- to the poly- 
disperse case, we plot in Fig. |9]the coexistence packing 
fraction of the cloud-nematic and shadow-smectic curves 
as a function of polydispersity for A = 1 and 2. The 
effect is weak when the temperature is below the tricriti- 
cal point of the monodisperse case, since the density gap 
hardly changes with A. By contrast, for temperatures 
above the tricritical-point temperature of the monodis- 
perse fluid, and as the polydispersity is increased from 
A = 0, the density gap opens up quite rapidly, mean- 
ing that the transition changes from first to second order 
even for a low degree of polydispersity. This is shown in 
Figs. M.c) and (d). 

We have also studied the effect of microsegregation, 
i.e of the spatial redistribution of particles with different 
sizes occurring along one period in the smectic phase in 
the coexisting shadow-smectic phase. For this purpose, 
the moment distribution profiles ma{z)^ defined previ- 
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FIG. 9: Packing fraction vs. polydispersity for attractive 
Yukawa fluid with: (a) and (c) A = 1; (b) and (d) A = 2. 
Temperatures are: (a) T* — 1.52, (b) 0.63 (both bellow corre- 
sponding tricritical points), (c) 2.50 and (d) 1.00 (both above 
corresponding tricritical points). 



ously, are useful. These functions are plotted in Fig. 
[TOTa). together with the distribution function h{z,r) as 
a function of r for fixed z in (b), and as a function of z 
for fixed r in (c), all for A = 0.294 and T* = 0.6. The 
radius distribution functions p{r) in the N (cloud) and 
S (shadow) phases are shown in Fig. [TOl d). The func- 
tions ma(z) reflect the fact that platelets are arranged in 
layers, with larger-sized platelets lying exactly at the lay- 
ers and with platelets being progressively smaller as one 
moves into the interstitial region. This microsegregation 
is more clearly appreciated in Fig. llOf b). where the size 
distribution function h{z,r) is plotted for selected val- 
ues of the z coordinate between the location of the layer 
z = and the intermediate distance z = d/2. At the 
layers, the maximum of the distribution occurs for larger 
values of the radius, whereas at the interstitial the maxi- 
mum is located at a rather low radius. The width of the 
distribution decreases from the layer to the interstitial. 
The shift in average size as one moves along a smectic 
period is again visible in Fig. [TUTc): of the three sizes 
considered, larger particles (r = 1.2) peak at the layers, 
and smaller ones at the interstitials, with intermediate 
sizes in intermediate positions. 

The size distribution functions of nematic and smectic 
phases are quite similar, with a slightly larger composi- 
tion of the larger-sized platelets in the smectic phase, but 
with almost identical widths. This trend is less and less 
apparent as the polydispersity A is increased: nematic 
and smectic distributions become more and more simi- 
lar, with mean radii that shift to lower values, widths 
that become broader and tails that decay more slowly. 
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FIG. 10: Radius distribution functions in the (shadow) coex- 
isting smectic phase for the attractive polydisperse Yukawa 
fluid with A = 2 at scaled temperature T* — 0.6 and de- 
gree of polydispersity A — 0.294. (a) Moments ma{z) for 
a = (continuous curve), 1 (dashed curve) and 2 (dotted 
curve), (b) Size distribution function h{z,r) as a function 
of r for different locations z along one smectic period; from 
right to left z/d — 0, 0.12, 0.24 and 0.5 (the smectic spac- 
ing is d = 1.194Lo). (c) Size distribution function h(z,r) as 
a function of z for different values of platelet radius; r = 1.2 
(continuous curve), 0.9 (dashed curve) and 0.5 (dotted curve), 
(d) Radius distribution functions p{r) for the coexisting ne- 
matic (at the cloud curve) and smectic (at the shadow curve). 



V. CONCLUSIONS 

In this work we have discussed the phase behaviour of 
suspensions made of platelets. We were motivated by the 
recent finding of smectic ordering in an aqueous suspen- 
sion of a— ZrP-based platelets of the same thickness but 
very polydisperse in diameter. The work has focused on 
the layered arrangement (smectic phase) of platelets aris- 
ing from a suspension of positionally disordered, but per- 
fectly orientationally ordered, particles (nematic phase). 
The high polydispersity in diameter prevents the forma- 
tion of the columnar phase, which has not been consid- 
ered in the present work. The assumption of parallel, 
perfectly oriented platelets should be valid when the as- 
pect (diameter-to-thickness) ratio is very high, which is 
the case in the a— ZrP samples and in many other discotic 
systems made from mineral materials. This assumption 
allows to use the powerful fundamental-measure density- 
functional theory for hard platelets, which is expected to 
accurately describe ordered arrangement of platelets, due 
to the emphasis of the theory on spatial correlations. Our 
first aim was to obtain a picture of how the degree of di- 
ameter polydispersity affects the nematic-smectic phase 
transition. Polydispersity tends to stabilise the smec- 
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FIG. 11: Evolution of the radius distribution function p(r) 
with increasing A for the nematic in the cloud curve (dashed 
curves) and the smectic in the shadow curve (continuous 
curves) at a scaled temperature T* — 0.625 and with an in- 
verse length A = 2. (a) A = 0.294, (b) 0.422 and (c) 0.523. 



tic phase, due to a microsegregation effect which allows 
larger and similar particles to populate the layers and 
pack more efficiently. The platelet volume fraction of the 
sample at the nematic-smectic transition can be reduced 
by 5% for a degree of diameter polydispersity of 30%, 
quite typical of the experimental system. The transition 
is found to be always continuous, regardless of the degree 
of polydispersity. 

Our second aim was to explore the consequences of soft 
interactions on the nematic-smectic transition. Suspen- 
sions of platelets can be obtained using different tech- 
niques, and both attractive and repulsive particles can 
be designed. We have contemplated both cases using a 
simple monotonic Yukawa-like interaction potential func- 
tion which can be positive or negative. The sign of the 



soft interactions profoundly affects the phase behaviour. 
For example, when platelets are identical (zero polydis- 
persity), soft repulsive forces always induce a continuous 
nematic-smectic transition, with the smectic phase sta- 
bilising with respect to the nematic for a decreasing range 
of the force. By contrast, when the interactions arc at- 
tractive, the transition becomes of the first order for tem- 
peratures below a tricritical temperature. In addition, 
when the interaction range is sufficiently long, there ap- 
pears a nematic demixing phenomenon, whereby two ne- 
matic phases of different platelet concentration appear in 
the suspension. Both nematic phases may coexist with a 
smectic phase at a triple-point temperature. In this case 
the smectic phase is slightly destabilised as the range of 
the force is decreased. Polydispersity changes this sce- 
nario to some extent. While nematic demixing is not 
much affected (at least for polydispersities < 30%), the 
tricritical temperature moves to very high values, and the 
nematic-smectic transition becomes of first order even for 
a small degree of polydispersity. 

The application of the present results to real experi- 
mental samples may be done after cautiously considering 
the type of interactions present in the platelets. Surface 
charges, counterions, hydration effects, etc. may play a 
role in the total balance between attractive and repulsive 
forces. Also, the interactions may have different sign with 
respect to distance. Another aspect, not contemplated 
in the present work, is the role of the isotropic phase, 
which is stable at lower packing fractions. The analy- 
sis of the isotropic phase requires a different theoretical 
framework in order to account for orientational degrees of 
freedom. For example, the fundamental-measure density 
functional for freely-rotating hard platelets with vanish- 
ing thickness [11] (adequately extended to the polydis- 
perse case) could be a good candidate for a reference 
system in perturbation theories for soft platelets. A 
complete picture would be obtained by incorporating the 
columnar phase, which could in principle be done within 
the present theory, but at a higher computational cost. 
The computation of the full phase diagram, and the study 
of how polydispersity in diameter and/or thickness, soft 
direct and solvent- induced interactions, etc., change the 
phase behaviour of platelet suspensions, is a big challenge 
from a theoretical point of view. This work represents a 
step toward that goal. 
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